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Abstract 

We compute the three-loop QCD corrections to the quark chromomagnetic mo- 
ment and thus obtain the matching coefficient and the anomalous dimension of the 
chromomagnetic interaction in HQET. As a byproduct we obtain the three-loop 
corrections to the quark anomalous magnetic moment. 
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1 Introduction 



We consider Quantum Chromodynamics (QCD) with ni light flavours and one 
heavy flavour Q. The interaction of a single heavy quark having momentum 
rnqv + k [mq is the on-shell mass and f ^ = 1) with gluons and light quarks 
in the situation when the residual momentum k <C mq (and momenta of 
light fields are also small) is described by the Heavy Quark Effective Theory 
(HQET) Lagrangian [HE] 



L = Q^iv ■ DQ., + {Ok + CcMOcM) + o(\ 
2mq \mf^^ 
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where 'f)Qv = Qv is the HQET quark field and D denotes the covariant deriva- 
tive (see the books [3.4J for more details). The kinetic energy operator 

Ok = -Q.DlQ, , D^ = D-v{vD), (2) 

does not renormalize and its coefficient is equal to 1, to all orders, due to 
reparametrization invariance [5]. The chromomagnetic interaction operator is 
defined as 

0cm '^QvO fii,(J^ Qv , (3) 

where G^i/ = QsG'^^j^ta is the gluon field strength tensor, ta is a colour matrix, 
as = gl/iyAn) the strong coupling and a^^ = |[7'^,7'^]- It is responsible for 
the violation of the heavy-quark spin symmetry and thus, e.g., for the B-B* 
mass splitting. The coefficient CcmifJ') is obtained by matching the scattering 
amplitudes of an on-shell heavy quark in an external chromomagnetic field, 
expanded in the momentum transfer q up to the linear term, in the full theory 
(QCD) and the effective theory (HQET). If all ffavours except Q are massless, 
all loop corrections in HQET vanish. It is most convenient to calculate the 
QCD scattering amplitude using the background ffeld method [6]. 

The chromomagnetic interaction coefficient CcmifJ') has been calculated at 
one-loop order in Ref . [1] . The one- loop anomalous dimension of the chro- 
momagnetic operator (JSj) follows from its fi dependence; it has also been found 
in Ref. [2]. The two-loop anomalous dimension has been obtained in Refs. [7|l8] . 
and the two- loop matching coefficient Ccm(/^) in Ref. pj. All orders of pertur- 
bation theory in the large-/5o limit were summed in Ref. [9]. The effect of 
a non-zero charm quark mass rric on the bottom-quark chromomagnetic in- 
teraction at two loops has been investigated in Ref. [TU]. In this paper (cf. 
Section [3]) we calculate 'jcm, as well as at three loops, provided that 

all light ffavours are massless. Using these results, we obtain the next-to-next- 
to-leading perturbative correction to the ratio 

R = f , 4 

which we discuss in Section m 

The mass difference between the vector and pseudo-scalar B meson is also 
often studied with lattice gauge theory simulations where non-perturbative 
results for the operator matrix element can be obtained (see, e.g., Ref. [11] for 
an introduction to lattice HQET and Ref. fT2] for a recent study). In principle 
there are two possibilities to make contact with the experimentally measured 
result: the matching can be performed perturbatively and non-pertubatively. 
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Fig. 1. Sample diagrams contributing to the quark chromomagnetic moment. Solid, 
curly and dotted lines denote quarks, gluons and ghosts, respectively. (8) represents 
the coupling of the background field. In the closed quark loops all flavours have to 
be considered. 

In the first case the perturbatively computed n-loop matching coefficient and 
the {n + l)-loop result for the corresponding anomalous dimension have to be 
combined with the non-perturbative lattice results. Currently this is done for 
n = 1 which according to Ref. |T2] induces an uncertainty of about 4%. With 
the new results of this paper this uncertainty can be significantly reduced. 
This topic is dicussed in Section |5l 

We also investigate the heavy-quark magnetic moments. They have not yet 
been measured experimentally, however, for the bottom and the lighter quarks 
there are upper limits from LEPl data [13]. For the bottom quark this limit 
is close to the Standard Model (SM) prediction including the two-loop QCD 
correction [T3]. Thus, a more precise measurement at a future linear collider 
should be able to determine the bottom-quark anomalous magnetic moment 
and probe possible deviations from the SM. 

The top-quark magnetic moment has not been measured so far. However, 
such a measurement would be very interesting since the top-quark couplings 
to photons or Z bosons are very sensitive to contributions from physics beyond 
the SM. Having this in mind, it is mandatory to have precise SM predictions 
for these couplings. In Section [6] we provide the three-loop QCD corrections 
to the coupling of the photon to heavy quarks. 



2 The Calculation 

To calculate the chromomagnetic moment we have to consider the quark- 
anti-quark-gluon vertex in the background-field formalism in QCD. Sample 
diagrams are depicted in Fig. [H When both the quark and anti-quark are 
on the (renormalized) mass shell and have physical polarizations, the vertex 
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= THa can be decomposed into two form factors, 

r^ = 7'^Fi(g2)--^a'^V2(g'), (5) 
2mQ 



where q = Pi — P2 '^^ the gluon momentum and pi and p2 are the momenta of 
the quark and anti-quark, respectively. 

The anomalous chromomagnetic moment is given by = Z^^F2{0), where 
is the quark wave function renormalization constant in the on-shell 
scheme. The total quark colour charge is given by -Z'^^Fi(O) = 1. Thus, -Fi(O) 
is the inverse of the on-shell wave function renormalization constant, which has 
been calculated to three-loops in Ref. [15] (see also Ref. [IS]). Therefore, the 
calculation of -Fi(O) provides a strong check on the correctness of our result. 

In order to extract the form factors, we use projection operators. They are 
conveniently obtained by introducing the momentum p = {pi + P2)/2, since 
p ■ q = 0. With this definition we have 



F.{q' 



2{d - 2){q^ - Ami) 



X Tr |(A + mg) (7. + ^^^T^P') + niQ) T^j , (6) 
2m^ 



f / + (d- 2)a'^ \ 1 

X TV |(A + ,n,) (7,. + ,4,.'_44) ] <'*^ + n 

Since the projector for F2 develops a pole for = 0, we cannot set q^ = from 
the beginning. Instead, we expand in q and keep all terms which are at most 
quadratic in q. In the final result the limit = can be taken. Due to the 
expansion in q all occurring integrals are on-shell propagator-type integrals. 

All Feynman diagrams are generated with QGRAF [TTj and the various topolo- 
gies are identified with the help of q2e and exp [TSfTQ] . In a next step the 
reduction of the various functions to so-called master integrals has to be 
achieved. For this step we use the so-called Laporta method [201I2T] which 
reduces the three-loop integrals to 19 master integrals. We use the implemen- 
tation of Laporta's algorithm in the program Crusher [22] . It is written in C++ 
and uses GiNaC [23] for simple manipulations like taking derivatives of polyno- 
mial quantities. In the practical implementation of the Laporta algorithm one 
of the most time-consuming operations is the simplification of the coefficients 
appearing in front of the individual integrals. This task is performed with the 
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help of Fermat [2l] where a special interface has been used (see Ref. |25j). 
The main features of the implementation are the automated generation of the 
integration-by-parts (IBP) identities [26j and a complete symmetrization of 
the diagrams. The master integrals are known from [IS] (see also comments in 
Ref. [ISj). To calculate the colour factors, we have used the program described 
in Ref. [2Z]Q 

The calculation is performed for an arbitrary gauge parameter in order to use 
its cancellation as a check. However, at three-loop level the expressions for the 
individual diagrams become very big. Thus we discard all terms with more 
than linear ^ dependence. This also concerns the factors 1/(1 — ^) appearing 
in the vertex of a background field with two quantum gluons, which does not 
appear in the usual formulation of QCD. If this is done, our final result is 
gauge-parameter independent up to terms which are quadratic in ^. Further- 
more, our calculation of Fi reproduces — including the gauge-dependent 
terms [T51ITB] . 



3 Chromomagnetic moment 



The renormalized scattering amplitude of an on-shell heavy quark with initial 
momentum pi = rngv and final momentum p2 = mqv — g in an external 
gluon field is given by the vertex ([5]) sandwiched between u{p2) and u{pi) and 
multiplied by 2'^^. We expand this amplitude in q up to linear terms, and 
re-express (relativistic) QCD spinors via HQET (non-relativistic) spinors: 



uirnqv + k) 



Uy{k) 



Then the QCD scattering amplitude reads 



2mc 



2mr 



a^^q,[l + ^i. 



taUv{0) 



(9) 



(we have used u{p2)^u{Pi) = 0). It must be reproduced by the HQET La- 
grangian of Eq. ([1]). If all flavours except Q are massless, all loop corrections 
vanish!^ and the scattering amplitude is given by the Born approximation: 



Uv{-q) 



2mQ 2mQ 



taM«(0) 



(10) 



^ We thank Philipp Kant for providing his interface for the program of Ref. [27] . 
^ We imply the use of Dimensional Regularization. 
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where = Z~J^{fi)CcmifJ') is the bare chromomagnetic interaction coef- 
ficient, and Zcmif^) is the MS renormahzation constant of the chromomag- 
netic operator given in Eq. ([3]). Note that both scattering amphtudes 
and fllUI) are renormahzed and hence ultraviolet-finite, however, both have in- 
frared divergences. These divergences are the same, because HQET has been 
constructed to reproduce the infrared behaviour of QCD. Vanishing loop cor- 
rection in HQET have ultraviolet and infrared divergences which cancel each 
other. The ultraviolet divergences of are removed by the infrared 

ones match those of 1 -|- /Xc (cf. Eq. (^). 



In order to find the chromomagnetic interaction coefficient Ccm{f^) in the 
HQET Lagrangian ([1]), we calculate the anomalous chromomagnetic moment 
He = -^2*^-^2(0) and re-express it in terms of a^^^\n), where the superscript 
denotes the number of active fiavours. Then 



where n/ = n/ — 1 is the number of light-quark fiavours, which are considered 
to be massless in our calculation, and rif is the total number of quark fiavours. 

The coupling constants in QCD (with ni + 1 fiavours) and a;^"''^(/i) 

in HQET (with ni fiavours) are related by [28] 



+ 



TT 



+ —L' 
V18 

a^)(/i)V 



36 



-TT 



+ 



n 



L+ (-L2 + — TT- )£ 



9 



77 



15 



31\ 



IL^ + _L + —7r' + — \Cp+ —L' --L + —tx' + 
4 8 48 32/ V12 9 144 



4 



43 
108 



\9 54 ' 



e + Oie'] \ + 0[a 



(12) 



where L = ln(/i^/mQ) and rriQ is the on-shell mass of the heavy quark. Cp = 
{Nl — l)/(2iVc) and Ca = N^. are the eigenvalues of the quadratic Casimir 
operators of the fundamental and adjoint representation for the SU(A^c) colour 
group, respectively. In the case of QCD we have Nc = 3 and Tp = 1/2. (n 
denotes Riemann's zeta function with integer argument n. 



The ultraviolet divergences contained in Zcm can be transformed into an 
anomalous dimension which is given by 
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7c 



din Zr 



a 




a 



in,)' 



TT 



17 

Ca I — Ca 

\36 



13 
36 



599 



65 \ 
216/ 



C3 + ^\C^. 



1728 
1 



1 JabcdJabcd 

2^ CfNf 



49 

2^^ ~ 96 



CaCf + —CaTftii 
36 



TfUi 



o JabcdJabcd 



-TT 



(13) 



where d"^^'^ 



and ci^'"^'^' 



are the symmetrized traces of four generators in the 



fundamental and adjoint representation, respectively (for SU(A^c), d'^'^'^d'^'^'^ 



{N^-l){Ny6N^ + 18)/{96N^), df^'^df'"^ = N,{Nyi){N^+6)/48). Nf = N, 
is the dimension of the fundamental representation. The two-loop result agrees 
with [71I5] . and the nf part of the three- loop one with [H]. 



Our result for Ccm reads 



Ccm(/^) = 1 + 



«("')(mQ) 



vr 



4 2 



Ca + -Cf 



a 



(ni) 



TT 



31\ 



-7r^ln2 + -C3 + — vr" ] C' 

2 4^ 12 16/ ^ 



,1. 1 2, 1 . 1 2 269 

+ -L + — TT^ In 2 - -Cs + — TT^ + 

* " 12 8^ 12 144 



CfCa 



+ H L H TT^ In 2 - -(3 H 

* 12 36 12 8 144 432 



^2 



25 



36 



13 



1 



299 



—CfTfUi + — - — L tt" C^T^n; 



24 



72 



36 



432 



+ 
+ 



^2 

- — 71 

3 



119 
~36 



CfTf + 



72 



216/ 



CaTf 



(14) 



The two-loop corrections were already calculated in Ref . [S] H] It is convenient 
to decompose the three-loop contribution in terms of the different colour struc- 
tures as 



4m ~ ^fffCf + XffaCfCa + XfaaCfCa + ^aaaCa 

^abcd^abcd 

+ ^dd AT [-^ffiCf + ^faiCfCa + ^aaiCa) TfUi 



^ Note that in Ref. [8j the term tt'^CaTf is not correct, since the O (e) term in the 
one-loop decoupling relation for (fT2|) has not been taken into account. 
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+ {XpFhCp + XfahCfCa + XaahCa) Tp 



+ {XfhhCf + XahhCa) n + {X^^m + . (15) 

Our results for the individual terms read 



20 5,4^ 5 2, 2^ 22 2, ^ 235, 103 

XpFF = —04 H In 2 TT^ In^ 2 tt^ In 2 U H vrVg 

3 18 18 3 24 72 ^'^ 

139 , 241, 23 2 101 

TT^ + Cs + — TT^ , 16 

2160 24 ^ 6 64 ' ^ ' 

/ 1 2, ^ 3 , 5 2 31\ , 35 35, 4^ 

XpFA^i --7rMn2H Cs H tt^ «4 ln^2 

V 8 16^ 48 64/ 3 72 

13 21 2 o 101 21 o 115^ 17 2. 209 4 847, 

+ — TT^ In^ 2 TT^ In 2 + U TrVg tt^ Cs 

9 24 6 ^9^ 2880 96 ^ 

9767 2 2803 

+ 1728^ -W ^^^^ 
1 .2 / 1 2i 1 . 1 2 337\ , 191 

191 4^ 169 21 2. 1745 2, ^ 265, 161 2. 

+ 432 - 216^ 2 + ^ - T^«= + ^» 

491 4 2951 17375 2 122971 

^ 10368^ ^ ^88"^^ ~ 3456 ^ ^ 10368 ' ^ ^ 

19,, 497 . /I 2, 1. 17 2 2917\ , 

432 1728 V48 32 ^'^ 576 3456/ 

8 1,4^ 1 2, 2. 317 2, ^ 925, 653 2. 

- -04 - - In^ 2 TT^ In^ 2 tt^ In 2 H Cs TrVg 

3 9 36 864 192^^ 864 

17 4 6079, 1585 2 1302797 , , 

TT (T-* H TT H fl9) 

810 1728^ 1296 186624' ^ ' 

^J/ = ^^^i^ - f «^ - ^ 2 + ^TT^ In^ 2 + ^TT^n 2 - IOC5 

97 2> 74 73, 151 2 5 

125 , , 

1 21 ^ 23 4 2, 23 2 88351 
TT^ In 2 TT^ + -C3 TT^ , (22) 

108 3240 3^ 72 10368' ^ ' 

35 „ 235 . /1> 1 2 715 \ , 4 1 , 4 

X^^i = + - -C3 + + L + -04 + — In^ 2 

"^^^ 864 864 V4^^ 144 1728/ 9 54 

1 21 2. 89 21 o 1 4 101. 35 2 

H TT^ In^ 2 TT^ In 2 + vr^ Cs + tt^ 

27 216 180 432^^ 864 
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236801 
46656 ' 



(23) 



X™ = lvr2 + — , (24) 

^" 27 324 ' ^ ' 

1 r3 13 ^2 1 r a 25 2 3535 ,^ , 

Xau = ^L -—L -^L + -C, + —rr (25) 

X™=-vr2- — , (26) 

108 162 ' ^ ' 



32 4,4^ 4 2, 2o 16 2, . 4 4 263^ 

—04 + - In^ 2 - -TT^ In^ 2 vr^ In 2 H Cs 

3 9 9 9 135 72 ^ 



11 2 2027 

H TT^H , (28) 

162 216 ^ ^ 

/ 1 2 119\ , 5,4 52,2 83 2, 

_?5c, + iA3-^.*-5H5Ic,-2??ll.^ + ?ilI, (29) 

24^^ 8 ^ 1440 864 19440 864 ' ^ ^ 

/ 5 2 149\ , 1,4 1 9, 2 1211 2, 

^AA/x = vr^ L + a4 + — In^ 2 tt^ In^ 2 + vr^ In 2 



TT _L + a4 H m 2: TT m 2: h 

288 864/ 24 24 432 

65 , 65 2 > 53 4 4423, 283429 2 

144^ 432 ^ 5184 1728 155520 
71965 



10368 ' 



(30) 



Xfhh--C3-j^^ (31) 

X^h,, = --Cs + — vr" + — , (32) 
Ahh 9'-^^ 270 972' ^ ^ 

= -Ivr^L + — TT^ - 3C3 - -TT^ + - , (33) 
3 270 27 3 ' ^ ^ 

= 16a4 + - In^ 2 - -tt^ In^ 2 - 2471^ In 2 + -(5 - -tt^Cs - — vr^ 
3 3 6^ 18 ^ 540 

4. 931 2 5 

with a4 = Li4(l/2). The nf part agrees with the result obtained in Ref. [9]. 

Substituting numerical values of the constants, we obtain, for the physical 
SU(3) colour group. 



CcU^q) = 1 + 0.6897a^)(mQ) + (2.2186 - 0.1938 n^) [a^)("^Q)] 

+ (11.079 - 1.7490 ri; + 0.0513 nf) (mg)] % O (af) 
= 1 + 0.6897 a("')(mQ) + (1.1626 /3o - 0.9786) Wj^'Xr 



[nic 
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+ (1.8468 + 0.3370 /?o - 3.813?) [a^;"\mQ)f + O (af) 



(35) 



with /?o = (11(^71/3 — 4Ti?n;/3)/4. The first-order l//3o result [9] contains the 
highest powers of /5o in each term. For rii = 4, for example, the coefficient of 



a^)(mQ) 



is 2.4221 - 0.9786 = 1.4435, and that of 
0.7020 - 3.8137 = 4.9039; the large-/3o approximation of Ref. [9 
includes the first terms in these sums, overestimates these two coefficients by 
68% and 63%, correspondingly. 



is 8.0155 + 
which only 



For the numerical evaluation of Ccmi^Q) and 7cm, we use the values rric = 
1.6 GeV, rrib = 4.7 GeV and rrit = 175 GeV. The number of light-quark 
flavours rii is three, four and five for the charm, bottom and top quark, 

defined with ni active flavours, from 



respectively. To evaluate a^J^^\mQ 



a 



(5), 



niz) = 0.118, we use the program RunDec j29] and obtain a^^^i 



rrir 



0.3348, 



nib 



0.2163 and a{^h 



0.1074, which leads to 



a 



cm \JT^c } 



Ccm{f^t) 



1 + 0.2309 + 0.1835 + 0.2362 = 1.6506 , 
1 + 0.1492 + 0.0676 + 0.0497 = 1.2664 , 
1 + 0.0741 + 0.0144 + 0.0045 = 1.0930 . 



(36) 
(37) 
(38) 



In the case of the charm quark, we see that the perturbative series does not 
converge which is probably connected to the relatively light scale of 1.6 GeV 
at which the strong coupling is evaluated. Thus, it seems that there are po- 
tentially large non-perturbative corrections. While the situation is better in 
the case of the bottom quark, the corrections are still very large. The three- 
loop correction amounts to about 30% of the one-loop contribution. For the 
top quark, we find that our new term contributes about 6% of the one-loop 
correction leading to a fairly reliable prediction of Ccm{fnt)- 

The numerical evaluation of the anomalous dimension ( fT3l) gives 



= 0.4775 a^"') + (0.4306 - 0.0549 rii) 
+ (0.8823 - 0.1472 ni - 0.0007 

For the individual quark flavours this leads to 



a 



(n,) 



a 



O a 



(39) 



7cm(mc) = 0.1599 + 0.0298 + 0.0163 = 0.2060 , 
lcm{mb) = 0.1033 + 0.0099 + 0.0029 = 0.1160 , 

Icmimt) = 0.0513 + 0.0018 + 0.0002 = 0.0533 . (40) 
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Thus, as far as the anomalous dimension is concerned the convergence be- 
haviour is acceptable even for the charm quark. 

These observations are in good agreement with the analysis of Ref. [9]. The 
chromomagnetic interaction coefficient Ccmifng) has the leading renormalon 
singularity (namely, a branch point) at the Borel parameter u = 1/2, quite 
close to the origin; it leads to a very fast growth of coefficients of the per- 
turbative series, ~ ?t,!(/5o/2)". It also means that the leading non-perturbative 
correction is only suppressed by the first power of I/uiq, and is thus impor- 
tant, especially for the charm quark. On the other hand, the perturbative 
series for the anomalous dimension has a finite radius of convergence. 



4 Application: heavy-meson mass splittings 



The most prominent physical effect caused by the chromomagnetic interaction 
is the mass splittings of hadronic doublets which are degenerate at mg = oo 
due to the heavy-quark spin symmetry. For example, for the bottom mesons 
B and B* one has [30], 

m% -ml = + O (^) , (41) 

where the index "(4)" means that we are considering rii = A fiavour HQET, 
and /ic'(4)(/u) is the matrix element of Ocm{lj) (cf- Eq. ([3])) over the ground- 
state meson. It is most natural to choose /i = mf, in Eq. (14T]) . because then Ccm 
contains no large logarithms. A similar formula can be written down for the D 
mesons where (/x) and /^^(s) (/i) appear in the corresponding expression for 
m\)* — mj^. The running of /^^•(^^(a^) governed by the anomalous dimension 
given in Eq. (fT3|) . Furthermore, it is necessary to relate the matrix elements 
in the two theories via the following decoupling relation 



/^G(4) ' 



rrir 



rrir. 



1 + Z2 



-Z3 



'a(^)(mc)' 



TT 



(42) 



with Z2 = — 71CaTp/432 [SHI]- We introduce the unknowrQ coefficient 
since it appears in estimates of higher order effects which are presented below. 

In the formulation with resummed logarithms one combines for consistency 
Eq. fH2l) with the two-loop result for Ccm and the three-loop anomalous di- 



Z3 can be calculated using 3-loop HQET integrals considered in |32|33j . 
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mension. This leads to the ratio R of Eq. (jl]) to the next-to-next-to-leading 
(NNL) order approximation. For later use we extend the formalism to NNNL 
order where we obtain 



R = x 



vr 



+ 

+ 



r2o + r2i(x2-l) + 



af {nib) ^ 



TT 



rso + r3i{x - l){x^ + X + 1) + -^{x - 1)^ + rir2o(x - 1) 

o 



+ rir2i(x - + 1) 



TT 



V ruhc 



(43) 



with 



X ■■ 



r2i 



gj^^ (^c) 

(4)/ \ 

as [mb) 



^^20 = c^^ini = 4) - 4^(ni = 3) + 2:2 



(2), 



-4^i(-. = 3) + ^^ 



+ 



7o 
4/3o 



+ Z2 



_72 ^ ATi ^ ^ _ /^A' 
7o /5o7o /5o V/^o, 



rso = cfi(nz = 4) - 43^(n, = 3) - c« (42i(n, = 4) - ^^^(nz = 3) + cis) 
+ 2:3, 

rsi = -c^^ini = 3) + c« (^^^(n, = 3) - ^2) - + z, + 



7o 
6/3o 



73 

7o 



/5o7o Polo \PoJ 7o /5o /5o /^o \f3o) 



(44) 



4^ denotes the coefficient of (^^(mQ )/7r)'^ in Ccm(/x = "^q)- The terms 2:2 
and stem from the decoupling of the matrix element and are introduced 
in Eq. (02]) and ^2 = [(2/9)Ca - {15/16)Cf]Tf stems from the decoupling 
of as (cf. Eq. (fT2l) ). 7„ are the coefficients of (as/Tr)""*"^ in the anomalous 
dimension and the coefficients of the f3 function are used in the form /?o = 
(11Ca/3 - 4TFn,/3)/4; see Refs. [51155] for the remaining A. 
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Both for 7„ and /3„ n/ = 4 active flavours have to be chosen. Let us mention 
that the next-to-leading (NL) order result of Eq. (143|) has been obtained in 
Ref. [7]. Inserting the numerical values given above and displaying the contri- 
butions from the individual orders separately, we find 



where the ellipses denote terms of higher order and power corrections. The 
term in square brackets is our estimate of the fourth order contribution, where 
we assume that the four-loop coefficient of the anomalous dimension is negligi- 
ble. This is justified by the rapid convergence of in the case of the bottom 
quark as can be seen in Eq. (HOl) . Furthermore, we set the unknown coefficient 
of Eq. (H2!) to zero which is a good approximation since it enters with a 
small coefficient. 

The experimental value is -Rexp = 0.88 with a negligible uncertainty. The 
NNL correction amounts to 10% of the LO contribution, however, it is larger 
than the NL one. Furthermore it is negative and thus increases the difference 
of the perturbative result and the experimental value. The estimated third- 
order correction is even larger than the NL and NNL one and contributes also 
with a negative sign. This indicates that the Aqcd / "^c correction may be quite 
substantial. 

It is interesting to consider the quantity R also without performing the resum- 
mation of the logarithms. In this way the three-loop result for the coefficient 
can be incorporated in a consistent way. The starting point is Eq. where 
quantities defined for rii = A are present in the numerator and the ones defined 
for n/ = 3 in the denominator. Using Eq. for the decoupling of the matrix 
element and running from = to /i = cancels /^c(4). Afterwards, we 
can replace af\'mc) by a^*^)(mb), using decoupling and renormahzation group 
running, and perform a consistent expansion of R in af\miy). As a result we 



R = 0.8517 - 0.0696 - 0.0908 + [-0.1285] + . . . 
= 0.6914+ [-0.1285] + ... , 



(45) 



obtain 



R = l-iol 



a 



vr 
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- 2/3ic« + (4/3o - 70)^2 - 72)/ + cfiim = 4) - 4^^(72, = 3) 



- 4^(4'i(^/ = 4) - c^l{ni = 3) + d^) + Z3 



+ a 



(46) 



where / = ln(m;,/mc). In our numerical evaluation we set the decoupling co- 
efficient ^3 to zero. Since similar three-loop decoupling effects are small we 
expect the same in our case. Furthermore, note that c^^ is numerically rather 
large. Inserting the numerical values yields 



R = l- 0.1113 - 0.0780 - 0.0755 + . . . = 0.7352 + 



(47) 



Comparing Eqs. (H5l) and (H71) . we find that the convergence of R without 
resummation behaves slightly better. However, the coefficients of the pertur- 
bative series are still large. 



5 Matching coefficient and renormalization group invariant quark 

mass 



In this section we would like to discuss the result of the matching coefficient 
in the form which is often used in lattice simulations of the mass difference 



m 



B* 



mg. In doing so we follow the procedure outlined in Ref. [37] . 



In lattice simulations one usually determines the renormalization group in- 
variant (RGI) matrix element of the operator Ocm which has to be multiplied 
by the corresponding matching coefficient. For its derivation one considers in 
a first step 



C„ 



2f3o exp 



TT 



7 



mag 



mag ^ 

7o 



2/3 2/3o 



a. 



(4J 



where m^, = mqijin^) is the scale invariant MS mass. The anomalous dimension 
7mag is given by a combination of 7cm («s), /?(««) and CcmifJ') and reads 



a 



a. 



s I mag I mag mag , 

7mag = — 70 + 7l — + 72 — 
TT V TT \ TT 



a. 



+ 



(49) 



:7cm + 2/3oc« (^)' + (4/?oc(^i + 2Ac« - 2/3o(cW)^) 



TT 
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The terms containing the /3-function stem from the running of a^. 

Since Cmag is still multiplied by the inverse pole mass, l/mq, in a further step 
the renormalization group invariant mass M^'^^ defined by 



M 



RGI 



TT 



7m. 
I30 



exp 



as{m.t) 



7m,o\ das 



In 



(50) 



can be used, where we introduced the quark mass anomalous dimension 7„ 



7m. I 7m, ~r 7m, 1 ~r 7m, 2 

TT \ TT 



+ 



(51) 



with 'jmfl = 'iCp/'^] the other coefficients can be found in Refs. j38]l39] . Using 
Eq. (l50l) in addition to the MS-on-shell relation the overall and logarithmic 
dependence on the pole mass can be replaced by the RGI mass. This procedure 
turns Cmag into Cgpin and one obtains an equation analog to Eq. (I48ll for Cgpin 
with the anomalous dimension given by 



7spin 7mag 



with 



(52) 



^mag 
^mag 
^mag 
^mag 



mag 

^7o 

mag 



:7r^-2/5ofci, 
: 7r' - Wiki + Po{2kl - 47oA;i - Ak2) , 
:73""" - W2ki + - ^loh - Ah) - l2Plhc'£ 



-Po (370^^1 -2kl - 671 /ei + Qkik2 - 670 /i;2 



6k' 



(53) 



The terms containing the /3-function stem from the transformation of rriQ 
to rfiQ and the coefficients ki are defined through mQ/fhQ{fhQ) = 1 + 
kias/n + k2a^/n'^ + ksa^/n^ + . . .. They can be found in Ref. [40j (see also 
Refs. piiii3l[T6] l 



In Fig. [2] the results for Cmag and Cspm 



are shown function of Aqcd/M^"^^ 
for ni = (left) and n/ = 4 (right) where the LO, NLO and NNLO results 
are shown. We also added an estimation for the NNNLO result by assuming a 
vanishing four-loop anomalous dimension which is motivated by the smallness 
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Fig. 2. Cmag and Cgpin as a function of Aqcd/'"^q for ri/ = (left) and = 4 (right). 
The upper group of hues corresponds to Cgpin while the lower shows Cmag- Inside 
these groups the dotted, dashed and solid lines show the LO, NLO and NNLO, 
respectively. The estimation for the NNNLO result is shown by the dash-dotted 
line. The region relevant for the bottom quark is indicated by the vertical line. 

of the higher order terms in Eq. ( |40|) . For the abscissa we choose Aqcd = 0.238 
which results from a lattice calculation for = [1^ and vary the value for 



For Kqcj:)/M^^^ 0.03, which is the range relevant for the bottom quark, one 
observes in the case of Cmag a relatively big shift when going from LO to NLO. 
However, the additional shifts after including the NNLO and the (estimated) 
NNNLO are smaller. The convergence improves significantly when going from 
Cmag to Cspin- In the case of the bottom quark the NLO corrections turn out 
to give a tiny contribution for rii = 0, however, also the NNLO and NNNLO 
results are quite small. Going to smaller quark masses one observes moderate 
corrections for Cgpin whereas for Cmag one has no convergence. We would like 
to mention that in the case n; = 4 there is basically no change in the behaviour 
of Cmag- However, as far as Cgpin is concerned one observes a moderate shift 
when including the NLO terms whereas the NNLO and NNNLO corrections 
are tiny. 



6 Magnetic moment 



The calculation of the anomalous magnetic moment of a heavy quark proceeds 
along the same lines as for the chromomagnetic moment. The only difference 
is that the external gluon is replaced by a photon. Our result reads 



Qq 



a 



(Uf) 



27r 



C. 



16 



(54) 

Qq is the charge of the heavy quark in terms of the positron charge. Note 
that the strong couphng in Eq. fl5^ is defined for Uf active flavours and it 
is evaluated at the scale fi = tuq. The two-loop contribution was already 
computed in Ref. [15]. Recently, it has also been obtained by considering the 
on-shell limit in the calculation of the off-shell form factor [1^ . We are in full 
agreement with Ref. [H] while we disagree with Ref. [45] by an overall factor 
of four in the coefficient of (as/vr)^. 

Our new three-loop term is given by 



+ 



+ 



a 



71 

317 1 



Cf 
1 



31 

16^ 12' 



-7r 



25 

+ ( ni + 

* 36 36 



144 8 4 g'-J / ^ 

119 1 



a 



TT 



''Q 



(3) /20 5,4^ 5 2,2^ 22 2, . 235, 103 2. 

a}^> = —04 + — In^ 2 TT^ In^ 2 vr^ In 2 Cs + TrVg 

^ V 3 18 18 3 24 72 ^'^ 

139 . 241, 23 2 lOr 



2160 24 ^^ 6 64 / 

20 5,4 49 2, 2 31 2, 185, 52. 

—04 + — In^ 2 vr^ In^ 2 + — vr^ In 2 (5 + —t^ Cz 

3 18 36 12 24 12 ^'^ 

35 4 113^ 1505 2 955\^2^ 

+ + Cs + ]CICa 

432 48 ^ 432 72 / ^ 

/5 5,4^ 11 2, 2o 25 2, . 65, 29 

+ -04 H In^ 2 TT^ In^ 2 H tt^ In 2 (5 H ti" Ca 

V3 72 18 8 32^^ 288 ^ 

103 4 1 463 2 ^31231 \ ^ 

H vr C'? TT H C-fG/1 

2880 2^ 216 2592 ) ^ ^ 

-a4 + - In^ 2 + -TT^ In^ 2 - -vr^ ln2 - — vr^ + 3(3 
3 9 9 3 216 ^ 



54 32 / ^ 



ia, + - In^ 2 + -Ti^ In^ 2 - -vr^ In 2 - — vr^ + -(3 
3 18 9 6 432 12^ 



77 2 2411\ 

H TT CpC^Tpn; 

216 324 / 

/32 4,4 4 2, 2 16 2, 4 4 263, 

162 864 y ^ 
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+ 



+ 



20 



32 



25 



—04 + — In^ 2 TT^ In^ 2 tt^ In 2 + — Cs - -vrYg 

3 18 18 9 24^^ 8 ^ 



143 



241 



2160 36 ^ 

1 2/ 317 
vr (n/ + 1) + T^ni 



-TT 



27 



4 



-Cs vr^ 

3^ 135 



324 
943 
324 



1375 

648 
61 

" 162 

JabcJabc 
^ Tn2 I "F "F V 
Upl p -\ yi-s^ 



CpTpni 



(55) 



where d'^'^ is the symmetrized trace of three generators in the fundamental 
representation (for SU(A^c), df^df" = {N^ - 1){N^ -4:)/{lQNc)). X,in denotes 
the contribution from singlet diagrams (cf. Fig. [T](g)). It is given by 



Xsin = 16a4 + - In^ 2 - -tt^ In^ 2 



3 

3^=^+ 54 



3 

931 2 



247r2ln2 + -C5 
6 



5 

18 



41 
540' 



-TT 



9 ' 



(56) 



where we only include the contribution from diagrams with closed heavy-quark 
loops. In principle there are contributions from diagrams with massless quarks 
as well. However, within perturbation theory they are divergent. Their eval- 
uation either requires non-perturbative methods or phenomenological models 
describing the interaction of light mesons in intermediate statesiU This is in 
analogy to the light-by-light contribution to the anomalous magnetic moment 
of the muon [46j, which contains logarithms of the electron mass. 



It is possible to obtain the known three-loop result for the electron anomalous 
magnetic moment from the expressions in Eqs. (!M|) - (!56|) by setting Cp = 
Tp = 1, Ca = 0, Np = 1, df'^df^ = 1 and = 0. This result was first 
obtained in analytical form in Ref. [20] and confirmed in Ref. [15]. 



Let us evaluate the quark magnetic moment numerically for charm, bottom 
and top quarks. Inserting the numerical values for the coefficients, we find 



^ = 0.2122 (mg) + (0.8417 - 0.0469 rii) [a^"/) (mg) 

Qq ^ 

+ (4.5763 - 0.5856 ni + 0.0145 nf) [^^"^^(mQ)] ^ + O (a^) . (57) 



Using a^^h 



rrir 



0.3378, ai^\mb) = 0.2169 and af\mt) = 0.1075 we obtain 



ac= 0.0478 + 0.0533 + 0.0758 = 0.1770, 

^ In principle it is possible to calculate the contributions from massive charm and 
bottom quarks to ah and a^, respectively. However, since the integrals involved 
contain two mass scales, this is beyond the scope of this work. 
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ttb = -0.0153 - 0.0103 - 0.0084 = -0.0340 , 
at= 0.0152 + 0.0047 + 0.0017 = 0.0215. 



(58) 



The pattern is very similar to the chromomagnetic moment of the heavy quark: 
no convergence is observed for the charm quark, the corrections for the bottom 
quark are large and amount to more than 50% of the one-loop contribution, 
and in the top quark case, we find that our new term gives a 10% contribu- 
tion which provides quite some confidence that the uncertainty of the final 
prediction for at is small. 

As already mentioned in the Introduction the LEPl bound of Ref. [13] for 
the bottom quark — at/Qb < 0.045 (68%C.L. ) — was found to be saturated 
by the two-loop correction. It is therefore interesting to see what happens if 
we include our three-loop term. For this purpose, we have to evaluate for 
/i = mz- We find 



Since the three-loop correction is almost as large as the two-loop one this 
overshoots the bound by about 25%. In this context we want to mention again 
that the contributions from closed light-quark loops could not be included in 
our calculation. These contributions might decrease the three-loop correction. 
In any case, a more precise measurement of would certainly be interesting. 



7 Conclusion 

Our main results are the three-loop anomalous dimension of the HQET 
chromomagnetic-interaction operator f[T5]l and the three-loop coefficient of this 
operator in the HQET Lagrangian f[T^ . They can be used for the investigation 
of various effects of the heavy-quark spin symmetry violation (e.g., mass split- 
tings) using continuum or lattice techniques. Furthermore, we have obtained 
the three-loop anomalous magnetic moments of heavy quarks fl54l) . This result 
does not include the light-quark light-by-light contribution which cannot be 
calculated perturbatively. 
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